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The presence of approximate electron number conservation and fj,-r permutation symmetry of 
5*2 is shown to naturally provide bilarge neutrino mixing. First, the bimaximal neutrino mixing 
together with [7 e 3=0 is guaranteed to appear owing to S2 and ,then, the bilarge neutrino mixing 
together with \ U e 3\ <0 arises as a result of tiny violation of S2. The observed mass hierarchy of ArriQ 
<C Arriatm is subject to another tiny violation of the electron number conservation. This scenario 
is realized in a specific model based on SU(3)l x U(1)n with two-loop radiative mechanism for 
neutrino masses. The radiative effects from heavy leptons contained in lepton triplets generate the 
bimaximal structure and those from charged leptons, which break S2, generate the bilarge structure 
together with |C/ e 3| -CI- To suppress dangerous flavor-changing neutral current interactions due to 
Higgs exchanges especially for quarks, this S2 symmetry is extended to a discrete Zg symmetry, 
which also ensures the absence of one-loop radiative mechanism. 

PACS numbers: 12.60.-i, 13.15.+g, 14.60.Pq, 14.60.St 



I. INTRODUCTION 

Recent observations of terrestrial neutrino oscillations at K2K and KamLAND 2] are proved to be consistent 
with the observations of atmospheric and solar neutrino oscillations 0, 3 ■ These oscillations are characterized by the 
squared mass differences and mixing angles, which are respectively given by Am 2 atm - 2 x 1CT 3 eV 2 and sin 2 28 atm ~ 1 
for atmospheric neutrinos and Ato| ~7x 1(T 5 eV 2 and sin 2 20q ~ 0.8 for solar neutrinos These mixings can be 
theoretically explained by the oscillations among three known neutrinos, v e ^ T , if they are massive [||. In fact, the 
experimental data indicate that the atmospheric and solar neutrino oscillations can, respectively, originate from the 
v^-v r mixing and the Ve-v^ mixing, which are both almost maximal. The mass hierarchy of Am 2 t 3> Attiq as well 
as the large mixing angles suggest that the neutrino mass matrix has almost bimaximal structure 0, II] ■ 

To accommodate massive neutrinos with tiny masses of C^lCr 1 ) eV (~ y/E mj fm eV), there are two theoretical 
ideas: one is the seesaw mechanism 0,0] and the other is the radiative mechanism 0,^3- While to understand 
the observed patterns of the neutrino mixings, let us consider, as the first approximation, the maximal neutrino 
mixings. To realize these maximal mixings, we often demand specific relations among neutrino masses, for example, 
TCiv^v lt — m ve"T ano Vor m v v =m v ^ Vr , where m ViVj (i,j=e,(j,,T) stands for the mass for Vi-Vj. To naturally ensure 
the presence of such specific relations, one may invoke a certain symmetry such as a U(l) symmetry based on 
L e — — L T (= L') [l3| for the maximal solar neutrino mixing or a fi-T permutation symmetry for the maximal 
atmospheric neutrino mixing 0, 0, 0, . 

In this article, we consider the neutrino oscillation based on the radiative mechanism. It has been argued that the 
original version of the Zee's one- loop radiative mechanism .11] failed to explain the favorable LMA solution for solar 
neutrinos, namely, the significant deviation from the maximal solar neutrino mixing |l8j . To get around the difficulty, 
one has to extend the original Zee's framework 0,^1 by 1) permitting the second Higgs scalar (introduced in the Zee 
model) to couple to the leptons, 2) employing a triplet Higgs scalar, 3) including a sterile neutrino or 4) introducing 
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two-loop radiative effects with the absence of one-loop radiative ones. We are concentrating on examining two-loop 
radiative effects [l2| since there are few theoretical models based on two-loop radiative mechanism that explain the 
observed properties of the current neutrino oscillations. The advantage of utilizing two-loop radiative effects lies in the 
fact that Majorana neutrino masses of O(10 _1 ) eV can be naturally generated without strong fine-tuning of various 
couplings. Denoting the lepton-number violating coupling by / inherent to Majorana neutrino masses, we roughly 
estimate the one-loop effect to be /m 2 ^/167r 2 M 2 and the two-loop effect to be / 2 m 2 /i/(167r 2 ) 2 M 2 for /i ~ v wea k 
(=(2V2G F )- 1 / 2 = 174 GeV) as the mass scale of the standard model and M ~ 1 TeV as the mass scale beyond the 
standard model. The estimated mass becomes ©(lCT 1 ) eV for / - 1CT 4 in the one-loop effect and / ~ 10 1 in the 
two-loop effect. The very small magnitude of / ~ 10 -4 for the one-loop case entails the strong fine-tuning while the 
magnitude of / ~ 0.1 for the two-loop case can be naturally of order of e, where e stands for the electric charge. In 
the present analysis, all the new couplings / beyond the standard-model ones are set to be ~ e. The pattern of the 
neutrino mixings is also naturally explained by the use of the fx-T permutation symmetry. 

To respect the power of the permutation symmetry (without fine-tunings of couplings to meet an "artificial" 
permutation symmetry) , since ordinary charged leptons badly break the permutation symmetry by their masses, new 
sources, whose interactions can preserve the permutation symmetry, are required. As the simplest extension of the 
standard model, we choose an SU(3)l X U(1)n model [2£j, which contains heavy leptons in lepton triplets as new 
sources. SU(3)l x U(1)n models are known to exhibit the intriguing aspect that 

1. it predicts three families of quarks and leptons if the anomaly free condition on SU(3)l X U(1)n and the 
asymptotic free condition on SU(3) C are imposed; 

2. the need of the three families of quarks and leptons group-theoretically arises because the anomalies from three 
families of leptons and one family of quarks that are SC/(3)z,-triplets are cancelled by those from two families 
of quarks that are iS , /7(3)i-antitriplets. 

Furthermore, these models naturally provide the radiative mechanism for neutrinos of the Majorana type because the 
standard Higgs doublet together with a Zee scalar specific to the one-loop radiative mechanism is regarded as the 
Higgs triplet |2l], |22, |23j . Our SU(3)l x U(1)n model includes three families of lepton triplets to be denoted b y 
— \v l L , t L , E l L ) T (i—e, fx, r), where (y % L , t L ) T is the standard doublet and E l L is negatively charged heavy leptons |24l. 
It can be argued that to activate a two- loop radiative mechanism requires a Higgs scalar of £ + , £ + ) T |23| . 

It will be shown that Ea. (|13fl is generated by leading two- loop radiative effects due to the heavy-lepton exchanges 
while the charged-lepton exchanges add less-dominant two-loop effects that give a slight deviation of a tm from the 
maximal value of 9 atm = tt/A and a significant contribution on Q to give sin 2 20q ~0.8. 

In the next section, we discuss the possible pattern of neutrino mass matrix, which is consistent with the observed 
neutrino oscillations. In Sec lIIII we describe the detail of the SU(3)l X U(1)n model including Yukawa and Higgs 
interactions as well as the symmetry structure summarized in Sec[n] In Sec lllll how to generate the neutrino masses 
and mixings is discussed and the observed properties of neutrino oscillations are shown to be well explained. The 
final section is devoted to summary. 



To see the appearance of underlying symmetries that approximately describe the observed neutrino oscillations, we 
show the parameterization of masses and mixing angles in terms of neutrino masses given by a matrix, M^, 



II. NEUTRINO MASS MATRIX 
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The digitalization of M v is performed by the mixing matrix of Umns B defined by 
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with 



The masses, 171^2,3, and mixing angles, #12,23,13, for 0q=0\ 2 and # a *m = #23, are computed to be 
mi = c 2 12 Xi + s 2 12 X 2 - 2c\ 2 s\ 2 A, m 2 = s 2 2 Ai + c 2 2 A 2 + 2C12S12A 

m 3 = sj 3 d + cj 3 f + 2c 23 s 23 e, (4) 
sin 2 2# 12 = ^^, sin 2 20 23 = 1 - 4A 2 ^c 2 3 = 1 + A, s 2 23 = i - 

, 0(3 S23^+C23C , . 

tan2#i 3 = 2 — (5) 

m 3 - Ai 

\ \ 2 j , 2 <r o C 2 3&— S23C . . 

Ai=o, A 2 = c 23 d+ s 23 / - 2c 23 S23e, A = , (6) 

Cl3 

x = V~2^-±, A = ^!M^, (7) 
A 4e- V2(6 + (tc)si 3 

for A 2 <Cl, where c a = cos# a and s a = sin6* a for a=12, 23, 13, and a — ±1 for S23 = =b j S23 1 - 

Since sin 2 #13 <Cl is reported, let us first set sin#i 3 =0. For the maximal atmospheric neutrino mixing char- 
acterized by A=0, the relation of d=f (corresponding to m v ^ v =m Vl .v T ) is required and can arise as a result of the 
fi-r permutation symmetry while the maximal solar neutrino mixing characterized by x=0 arises from a=d=e=f=0, 
indicating the U conservation. Our strategy used in this article to realize the bilarge neutrino mixing is to start with 
the maximal atmospheric neutrino mixing due to the /Lt-r permutation symmetry. This symmetry allows us to set 

a = 0, c = ±b, f = d, (8) 

which are subject to 

1. the L e conservation for a=0; 

2. a tiny breaking of the L e conservation for 6^0; 

3. the ji-T permutation symmetry for c=±b and f=d. 

It should be noted that a further possible relation of d = f = e (or d = / = — e) may arise from the Majorana mass 
term of [v T + v^) 2 (or (v T - ) 2 ) . 

Restricting ourselves to the case of d = / = e, thus utilizing (y T + v^) 2 , we find the following situations on the 
masses and mixing angles of neutrinos. In the case of C23 > and S23 > 0, the mass matrix with c = b yields 

mi = m 2 = 0, m 3 = 2d, sin 2 26> 23 = 1, tan26>i 3 = V2^, (9) 
where #12 is undetermined, while the mass matrix with c = — b yields 

|mi| = \m 2 \ = V2b, m 3 = 2d, sin 2 20 23 = sin 2 20 12 = 1, sin20i 3 = 0. (10) 
In the case of C23 > and S23 < 0, the mass matrix with c — b yields 

m 3 = 0, sin 2 20 23 = 1, sin 2 20 12 = — ^ (x = ^), sin2#i 3 = 0, (11) 

8 + x z b 

where mi and m 2 are given by Eq.Q with 0i 2 given here, while the mass matrix with c = —6 yields 

mi = sin 2 0i 2 d, m 2 = cos 2 12 d, m 3 = 0, sin 2 20 23 = 1, sin20i3 = 1, (12) 

where 0± 2 is undetermined. The gross feature of the observed neutrino oscillations is compatible with the zero-th 
order predictions in the cases of s 23 > with b = -c 0E! and of s 2 3 < with c = b . 
We choose the case of S23 > with c = — b that incorporates the bimaximal structure: 

/ b -b\ 

M v = b d d . (13) 
\-b d d J 

The hierarchy of \b\ <C \d\, which is the seed of Amg <C Am 2 tm realized after |mi| ^ \m 2 \, originates from the tiny 
violation of the electron number conservation. The observed property of sin 2 20 Q ~ 0.8 is explained by tiny violation 
of £2. In the next section, we realize this starting neutrino mass texture in a SU(3)l X U{1)n gauge model, where 
mass texture of heavy leptons as the third member of a lepton triplet plays an essential role. 
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III. MODEL 

The particles in our model are specified by quantum numbers placed as (SU(3)l, U(1)n). The hypercharge Y and 
the electric charge Q are given by Y—X s /^/3 + N and Q — (A 3 + Y)/2 respectively, where N/2 is the U(1)n number 
and A a are the Gell-Mann matrices with Tr(A°A b ) = 25 ab (a, b — 1,2, 8). Leptons are assigned to be: 

^l = {p\i\E i ) T L : (3,-2/3), r* r :(l,-l), E e R ^ + : (1,-1) , (14) 

where E R for j = (e, — ,+) are the mass eigenstates of the negatively charged heavy leptons and the superscripts, 
±, of E± R represent the chiral partners of the r ± p states of E±l to be defined by their Yukawa interactions while 
quarks are assigned to be: 

Ql = (d, -u, u') T L : (3*, 1/3) , Ql = (s, -c, c'f L : (3*, 1/3) , Ql = (t, b, b'f L : (3, 0) , 

u R ,c R ,t R : (1,2/3), d R ,s R ,b R : (1,-1/3), u' R , c' R : (1, 2/3) , b' R : (1,-1/3). (15) 

Another notations for Q l , Q' l =(d l ,u l ,u n ) L (i=e,p) and Q T =(u T ,d T ,d' T ) L , and for the corresponding right-handed 
quarks are also used. All gauge anomalies are cancelled by these quarks and leptons |24j. Higgs scalars are assigned 
to be: 

r,= (if\ V - ■■ (3,-2/3), p=(p+,p ,p°) T : (3,1/3), x = ( X + X°, X°f ■ (3, 1/3) , (16) 

which develop the following vacuum expectation values (VEV's): 

(0|r,|0) = K,0,0) T , (0|p|0) = (0, Wp ,0) T , (0|x|0) = (0,0,« X ) T , (17) 

where the orthogonal choice of these VEV's will be guaranteed by appropriate Higgs interactions introduced in Ea. (|26|) . 

Since the model respects the p,-r permutation symmetry of S2, in order to generate the phenomenologically con- 
sistent charged-lepton masses that do not respect the p-r permutation symmetry, the model employs the following 
^-antisymmetric Higgs scalars: 

p' = (p' + ,p'V°) T : (3,1/3), p" = (p"+,p" ,p"°) T : (3,1/3), (18) 
as well as SVsymmetric Higgs scalar, %': 

X' = (x' + ,X' ,X'°) T ■■ (3,1/3) (19) 

with 

<0|p'|0) = (0, Vl 0) T , (0|p"|0) = (0, V ',0) T , (0| X '|0)-(0,0,i; x T - (20) 

The two- loop radiative mechanism can be initiated by introducing three Higgs scalars denoted by k ++ and k' ++ 
as ^-symmetric states and £ as an ^-antisymmetric state: 

£= (C ++ ,C + ,H T : (3,4/3), ? = (e' ++ ,e' + ,f + ) T : (3,4/3),fc ++ ,fc'++ : (1,2), (21) 

where £' has L e =l and others have L e =0. 

Since our model contains quarks with the same charge, whose mass terms can be generated by p and \ between 
Q e L and down-type quarks and by p 1 and x 1 between Q^' T and up-type quarks, dangerous flavor-changing-neutral- 
currents (FCNC) interactions are generally induced at the phenomenologically unacceptable level j2^| • To avoid these 
interactions, Yukawa interactions must be constrained such that a triplet quark flavor gains a mass from only one 
Higgs field [23. The lepton sector also contains the similar FCNC problem because l l (i=e, p, r) and E l (i=e, ±) has 
the same charge. It is known that to impose such a constraint on FCNC is readily achieved by introducing a certain 
discrete symmetry. We combine the ^-permutation symmetry and the discrete symmetry for the FCNC-suppression 
into a Z% symmetry. Listed in Tableware the quantum numbers of S2, Zg, L and L e for leptons and Higgs scalars in 
our discussions, where ip±L = (ipl ± V'l)/'^, ^±-R = ( T R ^ Mfi)/ 1 /^ and £2 is shown for comparison. In this table, we 
have omitted the quark sector of the model since Zg of quarks can be easily adjusted to respect the Zg, conservation 
for the given Z% of Higgs scalars to reach quark interactions in Ea. (|22|l . 
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The Yukawa interactions for leptons are now described by Cy: 



Cy = W+L) C (fem +f(4>-LZ) 



1 



f^{i +R fl +R + -^ E {E. R fE_ 



k 



-+ 



f^(i. R f^ R + \f+ E {E +R fE +R 



,'++ 



with 



+ftW L pe R + ~¥ L {ftpt+R + fip't-R) + ~ L {gtp"t+R + 9jpi- R ) 
+f E ^l X E R + f+^xE +R + fsi^x'E-R 

+ ^ D R + P C U R + X c Ut) + 01 (VU R + pD T R + X D%) + (h.c.) (22) 



u r = E d r = E 4<4 u 'r = E ^ = (23) 



where /'s and g's denote the Yukawa couplings. The Higgs interactions are given by Hermitian terms composed of 
= V> Pi p' i p" i X, x'> £j ^ ++ ); which include the potential terms of V^ PX;PX 

V vpx = \ vp \ri x p\ 2 + X px \p x x| 2 + X xri \x x ?7| 2 , (24) 

y px = A^|pV|2 + a pp „|pVT + A xx <|xV| 2 (25) 

with the definition of (a x 6)" = e a ^ 'apb^ and by non-Hermitian terms in 

v = \\xu Wx)(?x') + tfrftfx)] + WfoW/O + \\w [faV)«V) + faV'XCV)] 

+P^vk ++ + m> vk '^vk' ++ + (h.c), (26) 

where fi's and A's denote mass scales and coupling constants, respectively. These Yukawa and Higgs interactions are 
invariant under 

• the .^-transformation, which is spontaneously broken; 

• the L e -transformation, which is not spontaneously broken but explicitly broken by ^rjk' ++ . 

The ^-permutation symmetry is preserved in the Yukawa interactions but explicitly broken by the Higgs interactions. 
The b- and c-terms in Eq.Q can be generated by interactions containing ^i]k' ++ . All other interactions are forbidden. 
We note that 

• the orthogonal choice of VEV's of T), p and x as in Ea. l|17|) is supported by V vpx if all A's are negative because 
V npx gets lowered if rj, p and x develop VEV's and one can choose VEV's such that (0|ryi|0) ^ 0, (0| p2 1 0) ^ 
and (0| X 3|0) ± 0; 

• the choice of VEV's of p', p" and x' as m Ea. (|20|l is supported by V px if all A's are negative; 

• a (ri'x)(£ ' p)-term, which is absent in V, would cause one-loop radiative effects such as in Fig^and similarly 
for the terms with p — > p' , p" and x X 1 ', 

• the suppression of FCNC for quarks is ensured by the absence of xD R and pD'^ for Q e L , of X*U R and p*U' R l for 
Q l L and of x^ R an d pE R for ip l L . 

Our model also respects 

• the L-conservation if A TO £ X ', A W £ P and A^ p '^ p " are absent; 

• the L '-conservation if p^nk' is also absent. 

Both the L- and L'-conservations are not spontaneously broken. 

Before discussing neutrino oscillations, we should mention another dangerous flavor-changing interactions of leptons 
due to the existence of p, p' and p" because the charged leptons can simultaneously couple to these Higgs scalars. 
These interactions are found to be suppressed down to the phenomenologically acceptable level because the suppression 
factor is at most m 2 . / v 2 ^^ from the incoming and outgoing vertices and provides enough suppression. Since the 
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approximate L e conservation is satisfied by our interactions, all L e -changing flavor interactions such as p — > ej are 
more suppressed. The detailed analyses have been described in Ref.01- The existence of heavy leptons and extra 
gauge bosons as well as heavy exotic quarks |2^| may also disturb well-established weak-interaction phenomenology. 
However, the additional contributions from these exotics will be well-suppressed because they are sufficiently heavy. 
Their masses controlled by the VEV's of x an d x' are taken to be ~ 2 TeV for the later analyses. 

The heavy lepton mass matrix is simply given by the diagonal masses computed to be mE<> = Jev x , m,E+ = f^ v x 
and rriE- = f7: v x'- On ^ ne other hand, the charged lepton mass matrix has the following non-diagonal form: 

/ mf \ 

Mt = mj" rof , (27) 

V vn\ T J 

where 

m T = feVp(= m e ), 

m T = \ {{ft + 97) v p ~ ft V - 9i v p"] > m T = \ {{ft - 97) v P + ffVpi - gjvp,,} , 

m T = \ {{ft - 97) v p - f7 v P' + 9ivp"} , m T e T = ~ [(/+ + gj) v p + f^Vp, + gjv p n] . (28) 

The diagonalized charged lepton matrix, Mf lag , are obtained by the use of two unitary matrices of U and V as 
shown in the Appendix. For the given mixing angles of a and (3 in Eq. HAlj) . the mass parameters in Ea. i|28|) can be 
determined by Ea. ljA5|l in terms of nip T . 



IV. NEUTRINO MASSES AND OSCILLATIONS 



The neutrino mass matrix is given by and M„ : 

M v = M* + M°, (29) 

where M„ (M„) arises from the two-loop radiative mechanism based on Fig|21(Fig[3J) for the heavy-lepton (charged- 
lepton) exchanges. The effective mass terms are expressed by 

(lp+LPX) 2 , i>eLPX ■ TP-LPX', ( 30 ) 

from Fig|2]and by 

(i>- L p + ^ +L p") X ■ ty-Lp + i>+Lp") x' + {X^ X'), 4>eLpX ■ (4>-LP + 4>+ L p") x! + (X «"> X'), (31) 

from Fig|3| where the product of the three particles reads abc=e a ^ 1 a a bpc 1 . After, Xi x\ P an( i p" acquire VEV's, 
where some of these contribute to masses of the heavy leptons and of the charged leptons, the heavy-lepton con- 
tributions supply main terms of neutrino masses, which are either ^-symmetric or ^-antisymmetric, while the 
charged-lepton contributions supply minor terms of neutrino masses, which are the admixtures of the ^-symmetric 
and SV antisymmetric contributions. These neutrino masses are expressed as: 

/ -r? r% \ ( Sml" 8m e J \ 

M* = -rf 1 1 J , = I §m e f Sm^ Sm£ T , (32) 

\ 1 1 J \ 8m e J Sm^ T Sml T J 

where and are calculated to be: 

= ~\fl (\p?p u p + \p>Hp"Vp>Vp») p^ k m E -I H (m|,m|-,m| ++ ) , 
t h _ .mf? IVyy ft" m 2 E+ Ih (rnl„m% + ,ml, ++ ) ^ 



ft P-ivk f k h m E - I H (m|,m|_,m^ ++ 
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with 



I H (m ,m E ,m k ) = 

l ( 2 2 2 2 2\ 

J {m 1 ,m 2 ,m 3 ,m 4 ;m k ) = 



J [m 2 E _ , m 2 ,m 2 E , m 2 _ ; — J fm|;, m 2 , 



m|_ , m 2 + ; to 2 



2 2 

to z — to; , 



G (to 2 , to|) G (m§, to; 



G(x,y) 



1 cc In (x/to^) - y In (y/m, 2 ) 



16tt 2 



x-y 



(34) 
(35) 
(36) 



under the assumption that the masses of k ++ and are heavy enough to neglect other masses in the integration. 
The masses used here are denoted by to? for the corresponding particle of <f>. The explicit form of J is given by 



J (a, b, c, d; e) 



d 4 k d 4 q 



1 



(2tt) 4 (2tt) 4 ( fc 2 _ a) (fc2 _ 6 j (g2 _ c) (g 2 _ d) fr k _ q f _ e ^ 
which can be evaluated via the following formula 



d 4 q 



(2tt) 4 (q 2 - a) (q 2 - b) (q 2 - c) ~ 16tt 2 



a In a 



blnb 



clnc 



(a - b) (a - c) {b - a) {b - c) (c - a) (c - b) 



(37) 



(38) 



To evaluate , it is convenient to express — Cy of Eq.J22J in terms of the neutrinos in weak interactions and the 
diagonalized charged leptons, ( e L, I^l, t l) T \ v=(y & j J , v ^l, v t l) t and £—(eji, [ir, tr) t , whose relevant terms are: 



v c U T f^Uip^ + + u c U T fyU^ l+ + e c V T f k V£k ++ + £ c V T f k ,Vek' ++ + ip e Mf a9 £ + (h.c), 



where 



ft 




1 



-1 -1 

^ 1 



fe 



f +e I 
4 = ^ Oil 
1 Oil 



fi 



+e ( 



0-10 

It is not difficult to reach an effective lagrangian corresponding to FigEl 



1 -1 
-1 1 



(39) 



(40) 



\xix' v x v x' 



E E ^a{U T ^U) mm ,{M f dlaa ) m , m , 

(<^,-R") = (£,fc),(£',fc') m,m' ,n,n'=e,fi,T 



where 



x (V T fKV) m , n ,I e (ml + ,m 2 K++ ) mini (M l diag ) n , nl (U T f,U) n , n v n , 

J(to 2 ,to 2 + ,to 2 ,to 2 _;to^J - J (to 2 ,to 2 + ,to 2 ,to| + ;to 2 -J 

h[m 0+ ,m K \ = g 2 

V v ) oh to _ — ml_L 



(41) 



(42) 



Ea. lfTTJl can be further calculated to be: 



-eff 



= -^fVt lX -nx(,x' v x v x'^kh (to| + ,to 2 ++ ) v c 



/0 



m 2 (c/3 + spf -to m to t (c 2 - sjA 
^ -m^mr (c| ~ s£) to 2 (cp - sp) 2 j 



16\/2 



x y 
x j (/. 
y 



(43) 



where the integral of Ig is defined as Ig (to 2 , to 2 ,) afc = 5 a fe/^ (to 2 , ti 2 ) with to 2 m t safely neglected in the denominators 
in the integral, and 



X = to 2 (c Q + s Q ) (cjg - s/?) 2 - to m to t (c a - s Q ) (c 2 , - s|) , 

y = TO M TO T (c Q + S a ) (c 2 - s|) - TO 2 (c a - S a ) (cp + Sp) 2 . 



(44) 
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Since 



£*ff = -\v c M°v, (45) 



the each entry of <5m„'s in Eg . (|32|) can be easily read off from Eg . 143(1 . 

To see the order of magnitudes of neutrino masses, we use the simplest choice of masses 



m 



1,h,t < m f r < m E-,E+ < m I+,e'+,€+,f'+( = m ^ < m l++,k'++- ( 46 ) 



This choice reduces the integrals of Ijj and Ig to 



I H {m 2 ,m%,m 2 k ) =-( — ) -^-^ In —fin — , £ ( m 2 , m j>) = - ( _ _^^ n _| m . (47) 



1 V 1 , ti? , to 2 



We also choose the following magnitudes of relevant parameters for neutrino masses: 

1. v p = Vpi = Vpn = v wea k/20 (~10 GeV) and v v = v wea k, where v p that controls the &-quark mass can be safely 
set to be ~10 GeV and v v that supplies the i-quark mass should be ~ v weak , since v v and v PtP > t p" are related to 
weak boson masses proportional to v 2 weak = £ aH v m ggs - 

2. v x — v' x — 10v wea k (~ 2 TeV) since \ and x' are the key fields for the symmetry breaking of SU(3)l xU(1)n — > 
SU(2) L x t/(l)y, leading to v XtX > > w luea fc- 

3. tyIjj- = v wea k, ?71e± = ev x (~ 0.5 TeV), mj- = m|- = m^i- = m^,- = v x (~ 2 TeV) and mk++,k'++ = v x /e (~ 6 
TeV) to satisfy Eq.gEJ. 

4. /i^jjfc = ev x with n^rjk 1 /jUf^fe ~ 0.1, which can be regarded as a natural relation |30| because the limit of ^'■qk' 
recovers the L e -conservation. 

5. /{^' = /j. = /fc/ = — Jfc = — fy = e and \p(,p,<qp' £p" ,vx£x' ~ e so that higher loop effects are safely neglected. 

This mass-setting provides sufficient suppression of exotic contributions in low-energy phenomenology. 

To see the magnitude of parameters for neutrino oscillations, we omit the terms proportional to (<C m T ) in 
Eq. (021 an d obtain that 

mf =3.3x 10- 2 - eV, rf - v^^l, (48) 

*m{J" = -6.6 x HT 3 ^ ( CjS + S/3 ) 2 eV, <5m^ = 6 - 6 x lO^ AW ( _ )2 ( + ^ (4Q) 
e ' V2 e 

where A = Aw^, = A r)(9 '£ P " and A^ = A^ x ^ x ' . By choosing the remaining parameters to be: 

- = 0.8, ^ = 2.0, ^2*1 = 0.12, (50) 
e A fj,£ n k 

so as to recover the observed data, we find using the estimates of Eas. 1(38)1 and l (35)l that, for the simplest case of 

Sa,f3 = 0, 



Ami = 2.2 x 10- 3 eV 2 , sin 2 29 atm = 0.96, 

Am% = 7.3 x 10~ 5 eV 2 , sin 2 26> = 0.79 (tan 2 6> = 0.37), 

1-2 



sin6»i 3 = 2.5 x 10~ 2 , (51) 
where cos 28 atm > (corresponding to A > 0) is chosen, and 

mi = 3.4 x 10~ 3 eV, m 2 = 9.2 x 10~ 3 eV, m 3 = 4.8 x 10~ 2 eV. (52) 
The mixing angles are determined by sin 26 atm = 1 - 4 A 2 and sin 2 26» = 8/(8 + x 2 ) with 

as in Eq©. The requirement of sin 2 6q ~ 0.8 (corresponding to x ~ v2) results in A2 ~ A. 

To see the dependence of the charged-lepton mixing angles, a and /3, in neutrino masses and mixings, we perform 
numerical estimation by keeping the terms proportional to in Eg. 1(43(1 and show the results in Fig|3|- Fig0 It can 
be stated that 
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1. As in Fig^I since the atmospheric neutrino mixing is dominated by the heavy- lepton contributions, the predic- 
tions of sin 2 29 atm and Am 2 tm are less-dependent on the charged lepton contributions, where their a-dependence 
hardly manifests itself in this graph because it reflects 1 < (c Q + s a ) < \[2 in Sniff 1 as in Eq. l|49|l . and similarly 
ffl3 is also less-dependent on the charged lepton contributions as in Fig[71 



2. As in FigO since the solar neutrino mixing heavily depends on the charged lepton contributions, the predictions 

4 



of sin 2 29q and Anig permitting the observed values vary with the magnitudes of and charged-lepton mixing 



angles; 

3. As in FigEl since the sin #13=0 from the heavy lepton contributions, its magnitude remain suppressed to be 
sin 013^0.03 induced by the charged lepton contributions. 

From Fig[3 we find that the effects of ^ in Ea. (|51|l lower (raise) the magnitude of sin 2 28 © (Am©) from 0.79 
(7.3 x 1(T 5 eV 2 ) to 0.76 (8.3 x 10~ 5 eV 2 ). We note that, in the simple-minded case of X e /X = 1 instead of X e /X = 2 
as in Eq. l|50f) . the neutrino oscillations are characterized by 

Am 2 tm = 2.3 x 10~ 3 eV 2 , sin 2 20 atm = 0.97, 

Aro| = 7.0 x 10~ 5 eV 2 , sin 2 20 Q = 0.85 (tan 2 Q = 0.44), 

sin 0i3 = (1.3- 1.4) x 10~ 2 , (54) 
for the mixing angles of the charged leptons at sin 2(3 = 0.75 with all values of sina(> 0) allowed. 



V. SUMMARY 

We have discussed how to theoretically understand the observed properties of neutrino oscillations, which can be 
summarized as: 

1. The atmospheric neutrino mixing is almost maximal and is characterized by sin 2 29 atm ~ 1; 

2. The solar neutrino mixing deviates from the maximal mixing and points to sin 2 20 Q ~ 0.8; 

3. The £/ e 3-mixing is suppressed; 

4. The squared mass differences show the hierarchy of Am 2 tm ^> Anig . 
Our answers consist of 

1. The atmospheric neutrino mixing is almost maximal because the dominant contributions to the neutrino mass 
matrix give 

/ b -b\ 

M„ = b d d , (55) 
\-b d d J 

as in Ea. (|13fl . which also give sin0i3 = 0, for C23 > and S23 > defined in Eq.J2J. 

2. The solar neutrino mixing deviates from the maximal mixing because less-dominant contributions perturb 
Eq.©, yielding sin 2 20 Q - 0.8; 

3. The J7 e 3-mixing is suppressed because the dominant contributions ensure the appearance of U e 3=0 and less- 
dominant contributions yield the suppressed U e s; 

4. The squared mass differences show the hierarchy of Am 2 tm ^> Am,Q because the tiny breaking of the electron 
number conservation allows the appearance of the suppressed Am 2 , . 

It should be noted that the following neutrino mass texture instead of Eq. (|55|l : 

/0 b b\ 

M' v =\ b d d , (56) 
\b d d j 

yields sin0i 3 (~ b[>/2d) 7^ 0. These observations are quite different from those based on the two-zero texture of 
neutrino masses |3l|. 

To realize the above scenario, we have adopted the model based on SU(3)l x U(1)n and have demonstrated that 
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1. The dominant contributions arise from the two- loop radiative effects of the heavy leptons; 

2. The less-dominant contributions arise from the two-loop radiative effects of the ordinary charged leptons. 

The model contains the specific Higgs scalars denoted by £ + , £ + ) T with L e =0 and by £' + ) T 

with L e = — 1, which initiate the two-loop radiative mechanism. These Higgs scalars can also initiate the one- loop 
radiative mechanism, which turns out to be forbidden by the discrete Zg symmetry that simultaneously forbids the 
appearance of dangerous flavor-changing interactions due to the direct Higgs exchanges inherent to this model. The 
dominant texture of Ea. (|55|l is subject to the ^-symmetric Yukawa interactions of the heavy leptons involving their 
diagonalized mass terms. The less-dominant charged lepton contributions receive the ^-breaking ones due to the 
existence of the ^-breaking charged-lepton masses. It is argued that, in Ea. l|55[l . the ^-contributions yield the d-term 
while the ^'-contributions yield the 6-term, which identically vanishes in the exact L e -conservation. 

The 5*2-permutation symmetry for /i-r can describe the gross feature of the observed atmospheric neutrino oscil- 
lations. Since the existing charged leptons obviously do not respect the ^-permutation symmetry, there should be 
other ingredients that respect the ^-permutation symmetry for the neutrino mixings. If we rely upon the seesaw 
mechanism instead of the radiative mechanism discussed in this article, it is expected that the type II seesaw mecha- 
nism | l(| is arranged to respect the ^-permutation symmetry while the ordinary seesaw mechanism 9] influenced by 
the charged leptons breaks the ^-conservation, which as in the present discussions provides the significant deviation 
of the solar neutrino mixing from the maximal one. 

APPENDIX A: CHARGED LEPTON MASSES 

In this appendix, the diagonalization of the lepton mass matrix Mi is demonstrated. The diagonal masses are 
obtained after the transformation by Mf lag = diag.(m e ,m^,m T ) = ujMgVi, where unitary matrices Ue and Vg are 
given by 

1 

V t = I eg sg \ (Al) 



with c a = cos a, etc., defined by 





( 1 


° 1 


u e = | 


c a 








c a ) 




+ (to7) 2 


- ml 


to 2 — to 2 




+ K T ) 2 


- ml 


TO 2 — TO 2 



-Sfi eg 



'-{my) 2 - 


K") 2 - 


h TO 2 


TO 2 


- ml 




-(to^) 2 - 


(to^ t ) 2 - 


h to 2 


TO 2 


- TO 2 





C0 = \ 5 5 , Sg = i/ - - , (A2) 

V m% — w m% — 

where to 2 + TO 2 = (TO^') 2 + (mJ r ) 2 +(m^ r ) 2 + (TOj M ) 2 is a obvious relation. The diagonal masses are computed to be: 

m 2 e = (to? 6 ) 2 , 

to 2 = i[(TOn 2 + (TOr) 2 +(m7) 2 + (TOn 2 -M 2 ], 

m 2 T = \ [(toD 2 + K" ) 2 + (to^) 2 + (mf ) 2 + M 2 ] , (A3) 

where 

M 4 = [{mYf - (toH 2 ] 2 + [{m7f ~ K T ) 2 ] 2 

+2(to[ t to^ t + mTmYf + 2(m T e T m T / + mfrof ) 2 . (A4) 

There are the following relations for non-diagonal and diagonal charged lepton masses: 

= S 2 m T + C m^, toJ t = C 2 m T + S' 2 to ai , 

[(c a s a C 2 - cgsgS 2 )m T - (cgsgC 2 - c a s a S 2 )m^\ , 



[{cgsgC 2 - c a s a S 2 )m T - (c a s a C 2 - CgSgS 2 )m^] , (A5) 
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where C 2 and S 2 are denned by 
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Table Captions 

TABLE [I] Particle contents with Sl2, and L e quantum numbers for quarks, leptons and Higgs scalars, where 
S2L = +{— ) denotes symmetric (antisymmetric) states and ui =1. 



Figure Captions 

FIG^ Forbidden one-lopp diagram. 

FIG[21 Heavy-lepton mediated two-loop diagrams. 

FIG |21 Charged- lepton mediated two- loop diagrams. 

FIG0] The dependences of sin 2 29 a t m (the black curves) and Am 2 tm (the gray curves) in functions of sin 2(3 for 
the given values of A^/A=1.0, 1.5, 2.0, where the dependence in sin a is so weak that it is covered by the line 
thickness. 

FIGE1 The dependences of sin 2 26*© (the black curves) and AtUq (the gray curves) in functions of sin 2(3 for the 
given values of sina=0 (thick curves), 1 (thin curves) and /A=1.0, 1.5, 2.0. 

FICHU The dependences of sin #13 as a function of sin 2(3 for the given values of sina=0 (thick curves), 1 (thin 
curves) and A£/A=1.0, 1.5, 2.0. 

FIG0 The dependences of mi. 2, 3 as a function of sin2/? for the given values of Xe/X—1.0, 1.5, 2.0, where the 
dependence in sin a is so weak that it is covered by the line thickness. 



TABLE I: Particle contents with Sl2, Zs and L e quantum numbers for quarks, leptons and Higgs scalars, where S2L = +(— ) 
denotes symmetric (antisymmetric) states and lo s =1. 
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FIG. 1: Forbidden one-lopp diagram. 
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FIG. 2: Heavy-lepton mediated two-loop diagrams. 
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FIG. 3: Charged-lepton mediated two-loop diagrams. 
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sin 2/5 



FIG. 4: The dependences of sin 2 26 a tm (the black curves) and Am 2 tm (the gray curves) in functions of sin 2/3 for the given 
values of Af/A=1.0, 1.5, 2.0, where the dependence in sin a is so weak that it is covered by the line thickness. 
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FIG. 5: The dependences of sin 2 28q (the black curves) and Ajtiq (the gray curves) in functions of sin 2/3 for the given values 
of sina=0 (thick curves), 1 (thin curves) and \i/X—1.0, 1.5, 2.0. 
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FIG. 6: The dependences of sin 613 as a function of sin 2/3 for the given values of sina=0 (thick curves), 1 (thin curves) and 
A*/A=1.0, 1.5,2.0. 
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sin 2/? 



FIG. 7: The dependences of 7711,2,3 as a function of sin 2(3 for the given values of A J >/A=1.0, 1.5, 2.0, where the dependence in 
sin a is so weak that it is covered by the line thickness. 



